Method and system for analyzing error probability of communication system using Gaussian Q function

ABSTRACT

Provided is a method and system for calculating a bivariate Gaussian Q function using a univariate Gaussian Q function to thereby analyze an error probability of a communication system. A method of calculating an error probability of a communication system may include: receiving encoded communication information from another communication system; restoring information desired to be transmitted by another user by decoding the encoded communication system; and calculating the error probability by analyzing the restored information using a bivariate Gaussian Q function, the error probability indicating a probability that the restored information is different from the information desired to be transmitted by the other user. The calculating may include calculating the bivariate Gaussian Q function using a univariate Gaussian Q function.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of Korean Patent Application No. 10-2009-0125956, filed on Dec. 17, 2009, and Korean Patent Application No. 10-2010-0029073, filed on Mar. 31, 2010, in the Korean Intellectual Property Office, the disclosures of which are incorporated herein by reference.

BACKGROUND

1. Field of the Invention

The present invention relates to a method and system for analyzing an error probability of a communication system by calculating a bivariate Gaussian Q function using a univariate Gaussian Q function.

2. Description of the Related Art

With developments of communication networks, there also increases a loss occurring when a communication system malfunctions due to an error.

Accordingly, to replace in advance a communication system having a relatively high error probability, an error probability of the communication system may be analyzed.

A bivariate Gaussian Q function may be used as a function to analyze the error probability of the communication system.

A conventional bivariate Gaussian Q function may be used to analyze the error probability of the communication system by employing a double integral structure as given by the following Equation 1:

$\begin{matrix} {{Q\left( {x,{y;\rho}} \right)} = {\frac{1}{2\pi\sqrt{1 - \rho^{2}}}{\int_{x}^{\infty}{\int_{y}^{\infty}{{\exp\left\lbrack {- \frac{u^{2} + v^{2} - {2\rho\;{uv}}}{2\left( {1 - \rho^{2}} \right)}} \right\rbrack}\ {\mathbb{d}u}\ {\mathbb{d}v}}}}}} & \left\lbrack {{Equation}\mspace{14mu} 1} \right\rbrack \end{matrix}$

However, an equation using the double integral structure may use a relatively large amount of calculations and thus, may require at least a performance for the above calculation. Accordingly, a system for analyzing the error probability of the communication system has been installed in an apparatus having a size greater than or equal to a predetermined level.

Accordingly, there is a desire for a method and system that enables even a mobile calculator of providing only simple calculations to calculate a bivariate Gaussian Q function without using a double integral.

SUMMARY

An aspect of the present invention provides a method and system that enables a low efficient calculator to analyze an error probability of a communication system by calculating a bivariate Gaussian Q function using a univariate Gaussian Q function that is a built-in function.

According to an aspect of the present invention, there is provided a method of calculating an error probability of a communication system, including: receiving encoded communication information from another communication system; restoring information desired to be transmitted by another user by decoding the encoded communication system; and calculating the error probability by analyzing the restored information using a bivariate Gaussian Q function, the error probability indicating a probability that the restored information is different from the information desired to be transmitted by the other user. The calculating may include calculating the bivariate Gaussian Q function using a univariate Gaussian Q function.

According to another aspect of the present invention, there is provided a method of calculating a bivariate Gaussian Q function, the method including: receiving values of x, y, and ρ that are variables of the bivariate Gaussian Q function; calculating values of ρ_(xy), ρ_(yx), and ρ that are variables of a bivariate Gaussian Q function; calculating the values of ρ_(xy), ρ_(yx), and δ_(xy) based on the values of x, y, and ρ; calculating a bivariate Gaussian function Q(x, 0; ρ_(xy)) using a univariate Gaussian Q function based on the values of x and ρ_(xy); calculating a bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function based on the values of y and ρ_(yx); and calculating the bivariate Gaussian Q function using Q(x, 0; ρ_(xy)), Q(y, 0; ρ_(yx)), and δ_(xy).

According to still another aspect of the present invention, there is provided a system for calculating an error probability of a communication system, including: a receiver to receive encoded communication information from another communication system; a decoder to restore information desired to be transmitted by another user by decoding the encoded communication system; and a bivariate Gaussian Q function analyzer to calculate the error probability by analyzing the restored information using a bivariate Gaussian Q function, the error probability indicating a probability that the restored information is different from the information desired to be transmitted by the other user. The bivariate Gaussian Q function analyzer may calculate the bivariate Gaussian Q function using a univariate Gaussian Q function.

According to yet another aspect of the present invention, there is provided a system for calculating a bivariate Gaussian Q function, the system including: a variable calculator to receive values of x, y, and ρ that are variables of the bivariate Gaussian Q function, and to calculate values of ρ_(xy), ρ_(yx), and δ_(xy) that are variables of a bivariate Gaussian Q function, based on the values of x, y, and ρ; a Q(x, 0; ρ_(xy)) calculator to calculate a bivariate Gaussian function Q(x, 0; ρ_(xy)) using a univariate Gaussian Q function based on the values of x and ρ_(xy); a Q(y, 0; ρ_(yx)) calculator to calculate a bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function based on the values of y and ρ_(yx); and a bivariate Gaussian Q function calculator to calculate the bivariate Gaussian Q function using Q(x, 0; ρ_(xy)), Q(y, 0; ρ_(yx)), and δ_(xy).

EFFECT

According to embodiments of the present invention, a bivariate Gaussian Q function may be calculated using a univariate Gaussian Q function that is a built-in function. Specifically, since the bivariate Gaussian Q function can be calculated without using a numerical integral, it is possible to calculate an error probability of a communication system using a low efficient calculator of providing only simple calculations.

BRIEF DESCRIPTION OF THE DRAWINGS

These and/or other aspects, features, and advantages of the invention will become apparent and more readily appreciated from the following description of exemplary embodiments, taken in conjunction with the accompanying drawings of which:

FIG. 1 is a block diagram illustrating a communication system using a Gaussian Q function according to an embodiment of the present invention;

FIG. 2 is a block diagram illustrating an example of a bivariate Gaussian Q function calculator of FIG. 1; and

FIG. 3 is a flowchart illustrating a method of calculating a bivariate Gaussian Q function according to an embodiment of the present invention.

DETAILED DESCRIPTION

Reference will now be made in detail to exemplary embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the like elements throughout. Exemplary embodiments are described below to explain the present invention by referring to the figures.

A bivariate Gaussian Q function that is in a double integral structure of above Equation 1 may be calculated using each of two different bivariate Gaussian Q functions as shown in the following Equation 2: Q(x,y;ρ)=Q(x,0,ρ_(xy))+Q(y,0,ρ_(yx))−δ_(xy)  [Equation 2]

Here, values of ρ_(xy), ρ_(yx), and δ_(xy) that are variables of a bivariate Gaussian Q functions of Equation 2 may be calculated by applying x, y, and ρ, which are variables of the bivariate Gaussian Q function, to the following Equation 3, Equation 4, and Equation 5:

$\begin{matrix} {\rho_{xy} = {{{sgn}(x)}{\left( {{\rho\; x} - y} \right)/\sqrt{x^{2} - {2\rho\;{xy}} + y^{2}}}}} & \left\lbrack {{Equation}\mspace{14mu} 3} \right\rbrack \\ {\rho_{yx} = {{{sgn}(y)}{\left( {{\rho\; y} - x} \right)/\sqrt{x^{2} - {2\rho\;{xy}} + y^{2}}}}} & \left\lbrack {{Equation}\mspace{14mu} 4} \right\rbrack \\ {\delta_{xy} = \left\{ \begin{matrix} {0;} & {{{xy} > {0\mspace{14mu}{or}\mspace{14mu}{xy}}} = {{{0\mspace{14mu}{and}\mspace{14mu} x} + y} \geq 0}} \\ {{1/2};} & {{AND}\mspace{14mu}{ELSE}} \end{matrix} \right.} & \left\lbrack {{Equation}\mspace{14mu} 5} \right\rbrack \end{matrix}$

Here, function sgn( ) of Equation 3 and Equation 4 corresponds to a sign function having a value of “1” when a value of ( ) is greater than or equal 0, and having a value of “−1” when a value of ( ) is less than 0.

Also, the bivariate Gaussian Q function may have a relationship as given by the following Equation 6: Q(−x,0;ρ)=½−Q(x,0;−ρ) Q(x,0;ρ)=Q(x)−Q(x,0;−ρ).  [Equation 6]

Here, the bivariate Gaussian function Q(x, 0; ρ) of Equation 6 may be calculated according to the following Equation 7:

$\begin{matrix} {{{{Q\left( {x,{0;\rho}} \right)} = {\int_{x}^{\infty}{\frac{\exp\left( {{- t^{2}}/2} \right)}{\sqrt{2\pi}}{Q\left( {- \frac{\rho\; t}{\sqrt{1 - \rho^{2}}}} \right)}\ {\mathbb{d}t}}}};}{{x \geq 0},{\rho < 0}}} & \left\lbrack {{Equation}\mspace{14mu} 7} \right\rbrack \end{matrix}$

When Equation 7 is substituted based on a univariate Gaussian Q function, it is possible to obtain an equation of calculating the bivariate Gaussian Q function using the univariate Gaussian Q function as given by the following Equation 8:

$\begin{matrix} {{{{Q\left( {x,{0;\rho}} \right)} \approx {{- \frac{\sqrt{1 - \rho^{2}}}{2\rho}}{\sum\limits_{i = 1}^{N}\;{\frac{a_{i}}{\sqrt{\beta_{i}}}{Q\left( {{- \frac{\rho\sqrt{\beta_{i}}}{\sqrt{1 - \rho^{2}}}}x} \right)}}}}};}{{x \geq 0},{\rho < 0}}} & \left\lbrack {{Equation}\mspace{14mu} 8} \right\rbrack \end{matrix}$

Here, a variable β_(i) of Equation 8 may satisfy the following Equation 9:

$\begin{matrix} {{{\beta_{i} = {\frac{1}{\sin^{2}\theta_{i}} + \frac{1 - \rho^{2}}{\rho^{2}}}};}{{{0 \leq \theta_{i}} = {{i\;{\pi/\left( {2N} \right)}} \leq {\pi/2}}},{i = 1},\ldots\mspace{14mu},N}} & \left\lbrack {{Equation}\mspace{14mu} 9} \right\rbrack \end{matrix}$

Specifically, a univariate Gaussian function Q(u) may be calculated using an error function erƒc(u) as given by the following Equation 10:

$\begin{matrix} {{Q(u)} = {\frac{1}{2}{{erfc}\left( {u/\sqrt{2}} \right)}}} & \left\lbrack {{Equation}\mspace{14mu} 10} \right\rbrack \end{matrix}$

The error function erƒc (u) of Equation 10 may be substituted with an equation of calculating an approximate value using an exponential function exp( ) as given by the following Equation 11:

$\begin{matrix} {{{{{erfc}(u)} \approx {\sum\limits_{i = 1}^{N}\;{a_{i}{\exp\left( {{- b_{i}}u^{2}} \right)}}}};}\left. N\rightarrow\infty \right.} & \left\lbrack {{Equation}\mspace{14mu} 11} \right\rbrack \end{matrix}$

Variables ai and bi of Equation 11 may satisfy the following Equation 12:

$\begin{matrix} {{{a_{i} = \frac{2\left( {\theta_{i} - \theta_{i - 1}} \right)}{\pi}},{b_{i} = \frac{1}{\sin^{2}\theta_{i}}}}{{{0 \leq \theta_{i}} = {{i\;{\pi/\left( {2N} \right)}} \leq {\pi/2}}},{i = 1},\ldots\mspace{14mu},{N.}}} & \left\lbrack {{Equation}\mspace{14mu} 12} \right\rbrack \end{matrix}$

FIG. 1 is a block diagram illustrating a communication system using a Gaussian Q function according to an embodiment of the present invention.

Referring to FIG. 1, the communication system using the Gaussian Q function may include an encoder 110, a transceiver 120, a decoder 130, and a bivariate Gaussian Q function analyzer 140.

The encoder 110 may encode information desired to be transmitted by a user using the communication system.

The transceiver 120 may transmit the encoded communication system to another communication system, or may receive encoded communication system from the other communication system.

The encoder 130 may restore information desired to be transmitted by another user by decoding the received communication system.

The bivariate Gaussian Q function analyzer 140 may calculate an error probability by analyzing the restored information using a bivariate Gaussian Q function. Here, the error probability indicates a probability that the restored information is different from the information desired to be transmitted by the other user.

As shown in FIG. 2, the bivariate Gaussian Q function analyzer 140 may include a variable calculator 210, a Q(x, 0; ρ_(xy)) calculator 220, a Q(y, 0; ρ_(yx)) calculator 230, and a bivariate Gaussian Q function calculator 240.

The variable calculator 210 may receive, from the decoder 130 of FIG. 1, values of x, y, and ρ that are variables of the bivariate Gaussian Q function, and may calculate values of ρ_(xy), ρ_(yx), and δ_(xy) that are variables of a bivariate Gaussian Q function, based on the values of x, y, and ρ. Here, the values of x, y, and ρ may correspond to parameter values of the restored information or communication information predetermined based on the communication system.

The variable calculator 210 may calculate the values of ρ_(xy), ρ_(yx), and δ_(xy) that are the variables of a bivariate Gaussian Q function by applying the received values of x, y, and ρ to Equation 3, Equation 4, and Equation 5, respectively.

The Q(x, 0; ρ_(xy)) calculator 220 may calculate a bivariate Gaussian function Q(x, 0; ρ_(xy)) using the univariate Gaussian Q function based on the values of x and ρ_(xy) calculated by the variable calculator 210.

When x is greater than or equal to 0 and ρ_(xy) calculated by the variable calculator 210 is less than 0, the Q(x, 0; ρ_(xy)) calculator 220 may calculate the bivariate Gaussian function Q(x, 0; ρ_(xy)) using the univariate Gaussian Q function according to Equation 8.

When x is greater than or equal to 0 and ρ_(xy) calculated by the variable calculator 210 is greater than or equal to 0, the Q(x, 0; ρ_(xy)) calculator 220 may calculate the bivariate Gaussian function Q(x, 0; ρ_(xy)) and the bivariate Gaussian function Q(−x, 0; ρ_(xy)) according to Equation 6.

The Q(y, 0; ρ_(yx)) calculator 230 may calculate a bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function based on the values of y and ρ_(yx) calculated by the variable calculator 210.

Specifically, when y is greater than or equal to 0, and ρ_(yx) is less than 0, the Q(y, 0; ρ_(yx)) calculator 230 may calculate the bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function according to Equation 8.

When y is greater than or equal to 0, and ρ_(yx) is greater than or equal to 0, the Q(y, 0; ρ_(yx)) calculator 230 may calculate the bivariate Gaussian function Q(y, 0; ρ_(yx)) and the bivariate Gaussian function Q(−y, 0; ρ_(yx)) according to Equation 6.

The bivariate Gaussian Q function calculator 240 may calculate the bivariate Gaussian Q function using Q(x, 0; ρ_(xy)) calculated by the Q(x, 0; ρ_(xy)) calculator 220, Q(y, 0; ρ_(yx)) calculated by the Q(y, 0; ρ_(yx)) calculator 230, and δ_(xy) calculated by the variable calculator 210.

Specifically, the bivariate Gaussian Q function calculator 240 may calculate, as the value of the bivariate Gaussian Q function, a value obtained by applying, to Equation 2, Q(x, 0; ρ_(xy)), Q(y, 0; ρ_(yx)), and δ_(xy), and by decreasing, by the value of δ_(xy), a sum of Q(x, 0; ρ_(xy)) and Q(y, 0; ρ_(yx)).

FIG. 3 is a flowchart illustrating a method of calculating a bivariate Gaussian Q function according to an embodiment of the present invention.

In operation S310, the variable calculator 210 may receive values of x, y, and ρ that are variables of a bivariate Gaussian Q function.

In operation S320, the variable calculator 210 may calculate values of ρ_(xy), ρ_(yx), and δ_(xy) that are variables of a bivariate Gaussian Q function, based on the values of x, y, and ρ.

Specifically, the variable calculator 210 may calculate ρ_(xy), ρ_(yx), and δ_(xy) that are the variables of a bivariate Gaussian Q function by applying the received values of x, y, and ρ to Equation 3, Equation 4, and Equation 5, respectively.

In operation S330, the Q(x, 0; ρ_(xy)) calculator 220 may verify whether x is greater than or equal to 0 and ρ_(xy) is less than 0.

In operation S340, when x is greater than or equal to 0 and ρ_(xy) is less than 0, the Q(x, 0; ρ_(xy)) calculator 220 may calculate the bivariate Gaussian function Q(x, 0; ρ_(xy)) using the univariate Gaussian Q function according to Equation 8.

In operation S350, when x is greater than or equal to 0 and ρ_(xy) calculated by the variable calculator 210 is greater than or equal to 0, the Q(x, 0; ρ_(xy)) calculator 220 may calculate the bivariate Gaussian function Q(x, 0; ρ_(xy)) and the bivariate Gaussian function Q(−x, 0; ρ_(xy)) according to Equation 6.

In operation S360, the Q(y, 0; ρ_(yx)) calculator 230 may verify whether y is greater than or equal to 0 and ρ_(yx) is less than 0.

In operation S370, when y is greater than or equal to 0, and ρ_(yx) is less than 0, the Q(y, 0; ρ_(yx)) calculator 230 may calculate the bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function according to Equation 8.

In operation S380, when y is greater than or equal to 0, and ρ_(yx) is greater than or equal to 0, the Q(y, 0; ρ_(yx)) calculator 230 may calculate the bivariate Gaussian function Q(y, 0; ρ_(yx)) and the bivariate Gaussian function Q(−y, 0; ρ_(yx)) according to Equation 6.

In operation S390, the bivariate Gaussian Q function calculator 240 may calculate the bivariate Gaussian Q function using Q(x, 0; ρ_(xy)) calculated in operation S340 or S350, Q(y, 0; ρ_(yx)) calculated in operation S370 or S380, and δ_(xy) calculated in operation S320.

Specifically, the bivariate Gaussian Q function calculator 240 may calculate, as the value of the bivariate Gaussian Q function, a value obtained by applying, to Equation 2, Q(x, 0; ρ_(xy)), Q(y, 0; ρ_(yx)), and δ_(xy), and by decreasing, by the value of δ_(xy), a sum of Q(x, 0; ρ_(xy)) and Q(y, 0; ρ_(yx)).

When a method of calculating a bivariate Gaussian Q function according to an embodiment of the present invention is applied to a conventional M-ary Phase Shift Keying (MPSK) symbol error probability equation P(E), P(E) may be modified as given by the following Equation 13: P(E)=1−Q[−√{square root over (2γ)} sin(π/M),−√{square root over (2γ)}sin(π/M);−cos(2π/M)]  [Equation 13]

Here, γ denotes a signal-to-noise ratio (SNR) per symbol and M denotes a magnitude of a symbol.

According to embodiments of the present invention, a bivariate Gaussian Q function calculation system may calculate a bivariate Gaussian Q function using a univariate Gaussian Q function that is a built-in function. Specifically, since the bivariate Gaussian Q function calculation system can calculate the bivariate Gaussian Q function without using a numerical integral, it is possible to calculate an error probability of a communication system using a mobile calculator or a low efficient calculator of providing only simple calculations.

Although a few exemplary embodiments of the present invention have been shown and described, the present invention is not limited to the described exemplary embodiments. Instead, it would be appreciated by those skilled in the art that changes may be made to these exemplary embodiments without departing from the principles and spirit of the invention, the scope of which is defined by the claims and their equivalents. 

1. A method of calculating an error probability of a communication system, comprising: receiving encoded communication information from another communication system; restoring information desired to be transmitted by another user by decoding the encoded communication information; and calculating the error probability by analyzing the restored information using a bivariate Gaussian Q function, the error probability indicating a probability that the restored information is different from the information desired to be transmitted by the other user, wherein the calculating comprises calculating the bivariate Gaussian Q function using a univariate Gaussian Q function.
 2. The method of claim 1, wherein the calculating comprises: receiving values of x, y, and ρ that are variables of the bivariate Gaussian Q function in the restored information; calculating values of ρ_(xy), ρ_(yx), and δ_(xy) that are variables of a bivariate Gaussian Q function; calculating the values of ρ_(xy), ρ_(yx), and δ_(xy) based on the values of x, y, and ρ; calculating a bivariate Gaussian function Q(x, 0; ρ_(xy)) using the univariate Gaussian Q function based on the values of x and ρ_(xy); calculating a bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function based on the values of y and ρ_(yx); and calculating the bivariate Gaussian Q function using Q(x, 0; ρ_(xy)), Q(y, 0; ρ_(yx)), and δ_(xy).
 3. The method of claim 2, wherein the calculating of Q(x, 0; ρ_(xy)) comprises: verifying whether x is greater than or equal to 0, and ρ_(xy) is less than 0; calculating the bivariate Gaussian function Q(x, 0; ρ_(xy)) using the univariate Gaussian Q function when x is greater than or equal to 0, and ρ_(xy) is less than 0; and calculating the bivariate Gaussian function Q(x, 0; ρ_(xy)) using the bivariate Gaussian Q function when x is greater than or equal to 0, and ρ_(xy) is greater than or equal to
 0. 4. The method of claim 3, wherein: the calculating of the bivariate Gaussian function Q(x, 0; ρ_(xy)) using the univariate Gaussian Q function comprises calculating the bivariate Gaussian function Q(x, 0; ρ_(xy)) using the univariate Gaussian Q function by employing ${{Q\left( {x,{0;\rho}} \right)} \approx {{- \frac{\sqrt{1 - \rho^{2}}}{2\rho}}{\sum\limits_{i = 1}^{N}\;{\frac{a_{i}}{\sqrt{\beta_{i}}}{Q\left( {{- \frac{\rho\sqrt{\beta_{i}}}{\sqrt{1 - \rho^{2}}}}x} \right)}}}}};$ x ≥ 0, ρ < 0, and variables α_(i), β_(i) , and θ_(i) are calculated by α_(i)=2(θ_(i)−θ_(i-1))/π; ${\beta_{i} = {\frac{1}{\sin^{2}\theta_{i}} + \frac{1 - \rho^{2}}{\rho^{2}}}};$ 0 ≤ θ_(i) = i π/(2N) ≤ π/2, i = 1, …  , N, N being a natural number.
 5. The method of claim 2, wherein the calculating of Q(y, 0; ρ_(yx)) comprises: verifying whether y is greater than or equal to 0, and ρ_(yx) is less than 0; calculating the bivariate Gaussian function Q(y,0; ρ_(yx)) using the univariate Gaussian Q function when y is greater than or equal to 0, and ρ_(yx) is less than 0; and calculating the bivariate Gaussian function Q(y, 0; ρ_(yx)) using the bivariate Gaussian Q function when y is greater than or equal to 0, and ρ_(yx) is greater than or equal to
 0. 6. The method of claim 5, wherein: the calculating of the bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function comprises calculating the bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function by employing ${{Q\left( {x,{0;\rho}} \right)} \approx {{- \frac{\sqrt{1 - \rho^{2}}}{2\rho}}{\sum\limits_{i = 1}^{N}\;{\frac{a_{i}}{\sqrt{\beta_{i}}}{Q\left( {{- \frac{\rho\sqrt{\beta_{i}}}{\sqrt{1 - \rho^{2}}}}x} \right)}}}}};$ x ≥ 0, ρ < 0, and variables α_(i), β_(i) and θ_(i) are calculated by α_(i)=2(θ_(i)−θ_(i-1))/π; ${\beta_{i} = {\frac{1}{\sin^{2}\theta_{i}} + \frac{1 - \rho^{2}}{\rho^{2}}}};$ 0 ≤ θ_(i) = i π/(2N) ≤ π/2, i = 1, …  , N, N being a natural number.
 7. The method of claim 2, wherein the calculating of the bivariate Gaussian Q function comprises calculating a value of the bivariate Gaussian Q function by decreasing a sum of Q(x, 0; ρ_(xy)) and Q(y, 0; ρ_(yx)) by the value of δ_(xy).
 8. A system for calculating an error probability of a communication system, comprising: a receiver to receive encoded communication information from another communication system; a decoder to restore information desired to be transmitted by another user by decoding the encoded communication information; and a bivariate Gaussian Q function analyzer to calculate the error probability by analyzing the restored information using a bivariate Gaussian Q function, the error probability indicating a probability that the restored information is different from the information desired to be transmitted by the other user, wherein the bivariate Gaussian Q function analyzer calculates the bivariate Gaussian Q function using a univariate Gaussian Q function.
 9. The system of claim 8, wherein the bivariate Gaussian Q function analyzer comprises: a variable calculator to receive values of x, y, and ρ that are variables of the bivariate Gaussian Q function, and to calculate values of ρ_(xy), ρ_(yx), and δ_(xy) that are variables of a bivariate Gaussian Q function, based on the values of x, y, and ρ; a Q(x, 0; ρ_(xy)) calculator to calculate a bivariate Gaussian function Q(x, 0; ρ_(xy)) using the univariate Gaussian Q function based on the values of x and ρ_(xy); a Q(y, 0; ρ_(yx)) calculator to calculate a bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function based on the values of y and ρ_(yx); and a bivariate Gaussian Q function calculator to calculate the bivariate Gaussian Q function using Q(x, 0; ρ_(xy)), Q(y, 0; ρ_(yx)), and δ_(xy).
 10. The system of claim 9, wherein the Q(x, 0; ρ_(xy)) calculator calculates the bivariate Gaussian function Q(x, 0; 92 _(xy)) using the univariate Gaussian Q function when x is greater than or equal to 0, and ρ_(xy) is less than 0, and calculates the bivariate Gaussian function Q(x, 0; ρ_(xy)) using the bivariate Gaussian Q function when x is greater than or equal to 0, and ρ_(xy) is greater than or equal to
 0. 11. The system of claim 10, wherein: the Q(x, 0; ρ_(xy)) calculator calculates the bivariate Gaussian function Q(x, 0; ρ_(xy)) using the univariate Gaussian Q function by employing ${{Q\left( {x,{0;\rho}} \right)} \approx {{- \frac{\sqrt{1 - \rho^{2}}}{2\rho}}{\sum\limits_{i = 1}^{N}\;{\frac{a_{i}}{\sqrt{\beta_{i}}}{Q\left( {{- \frac{\rho\sqrt{\beta_{i}}}{\sqrt{1 - \rho^{2}}}}x} \right)}}}}};$ x ≥ 0, ρ < 0, and variables α_(i), β_(i) and θ_(i) are calculated by α_(i)=2(θ_(i)−θ_(i-1))/π; ${\beta_{i} = {\frac{1}{\sin^{2}\theta_{i}} + \frac{1 - \rho^{2}}{\rho^{2}}}};$ 0 ≤ θ_(i) = i π/(2N) ≤ π/2, i = 1, …  , N, N being a natural number.
 12. The system of claim 9, wherein the Q(y, 0; ρ_(yx)) calculator calculates the bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function when y is greater than or equal to 0, and ρ_(yx) is less than 0, and calculates the bivariate Gaussian function Q(y, 0; ρ_(yx)) using the bivariate Gaussian Q function when y is greater than or equal to 0, and ρ_(yx) is greater than or equal to
 0. 13. The system of claim 12, wherein: the Q(y, 0; ρ_(yx)) calculator calculates the bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function by employing ${{Q\left( {x,{0;\rho}} \right)} \approx {{- \frac{\sqrt{1 - \rho^{2}}}{2\rho}}{\sum\limits_{i = 1}^{N}\;{\frac{a_{i}}{\sqrt{\beta_{i}}}{Q\left( {{- \frac{\rho\sqrt{\beta_{i}}}{\sqrt{1 - \rho^{2}}}}x} \right)}}}}};$ x ≥ 0, ρ < 0, and variables α_(i), β_(i) and θ_(i) are calculated by α_(i)=2(θ_(i)−θ_(i-1))/π; ${\beta_{i} = {\frac{1}{\sin^{2}\theta_{i}} + \frac{1 - \rho^{2}}{\rho^{2}}}};$ 0 ≤ θ_(i) = i π/(2N) ≤ π/2, i = 1, …  , N, N being a natural number.
 14. The system of claim 9, wherein the bivariate Gaussian Q function calculator calculates a value of the bivariate Gaussian Q function by decreasing a sum of Q(x, 0; ρ_(xy)) and Q(y, 0; ρ_(yx)) by the value of δ_(xy).
 15. The system of claim 9, wherein the values of x, y, and ρ that are the variables of the bivariate Gaussian Q function correspond to parameter values of the restored information or communication information predetermined based on the communication system.
 16. A method of calculating a bivariate Gaussian Q function for analyzing an error probability of a communication system, the method comprising: receiving values of x, y, and ρ from the communication system, x, y, and ρ being parameters of the communication system and being used as variables of the bivariate Gaussian Q function; calculating values of ρ_(xy), ρ_(yx), and δ_(xy) that are variables of a bivariate Gaussian Q function; calculating the values of ρ_(xy), ρ_(yx), and δ_(xy) based on the values of x, y, and ρ; calculating a bivariate Gaussian function Q(x, 0; ρ_(xy)) using a univariate Gaussian Q function based on the values of x and ρ_(xy); calculating a bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function based on the values of y and ρ_(yx); calculating the bivariate Gaussian Q function using Q(x, 0; ρ_(xy)), Q(y, 0; ρ_(yx)), and δ_(xy), and calculating the error probability for the communication system using the calculated bivariate Gaussian Q function.
 17. The method of claim 16, wherein the calculating of Q(x, 0; ρ_(xy)) comprises: verifying whether x is greater than or equal to 0, and p_(xy) is less than 0; calculating the bivariate Gaussian function Q(x, 0; ρ_(xy)) using the univariate Gaussian Q function when x is greater than or equal to 0, and ρ_(xy) is less than 0; and calculating the bivariate Gaussian function Q(x, 0; ρ_(xy)) using the bivariate Gaussian Q function when x is greater than or equal to 0, and ρ_(xy) is greater than or equal to
 0. 18. The method of claim 16, wherein the calculating of Q(y, 0; ρ_(yx)) comprises: verifying whether y is greater than or equal to 0, and ρ_(yx) is less than 0; calculating the bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function when y is greater than or equal to 0, and ρ_(yx) is less than 0; and calculating the bivariate Gaussian function Q(y, 0; 92 _(yx)) using the bivariate Gaussian Q function when y is greater than or equal to 0, and ρ_(yx) is greater than or equal to
 0. 19. A system for calculating a bivariate Gaussian Q function, the system comprising: a variable calculator to receive values of x, y, and ρ that are variables of the bivariate Gaussian Q function, and to calculate values of ρ_(xy), ρ_(yx), and δ_(xy) that are variables of a bivariate Gaussian Q function, based on the values of x, y, and ρ; a Q(x, 0; ρ_(xy)) calculator to calculate a bivariate Gaussian function Q(x, 0; ρ_(xy)) using a univariate Gaussian Q function based on the values of x and ρ_(xy); a Q(y, 0; ρ_(yx)) calculator to calculate a bivariate Gaussian function Q(y, 0; ρ_(yx)) using the univariate Gaussian Q function based on the values of y and ρ_(yx); and a bivariate Gaussian Q function calculator to calculate the bivariate Gaussian Q function using Q(x, 0; ρ_(xy)), Q(y, 0; ρ_(yx)), and δ_(xy). 